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ABSTRACT 

We report the first observational search for tau neutrinos (z/,-) from gamma- 
ray bursts (GRBs) using one of the Ashra light collectors. The earth-skimming 
Vj. technique of imaging Cherenkov r showers was applied as a detection method. 
We set stringent upper limits on the v^- fluence in PeV-EeV region for 3780 s 
(between 2.83 and 1.78 hr before) and another 3780 s (between 21.2 and 22.2 hr 
after) surrounding GRB081203A triggered by the Swift satellite. This first search 
for PeV-EeV Vr complements other experiments in energy range and method- 
ology, and suggests the prologue of "multi-particle astronomy" with a precise 
determination of time and location. 

Subject headings: gamma-ray burst: individual(GRB 081203A) — methods: ob- 
servational — neutrinos 



1. Introduction 

Gamma-ray bursts (GRBs) eject the most energetic outflows in the observed universe, 
with jets of material expanding relativistically into the surrounding interstellar matter with 
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a Lorentz factor V of 100 or more. Energy dissipation processes involving nonthermal in- 
teractions between particles are thought to play an important role in GRBs, but remain 
observationally unresolved. The detection of PeV-EeV neutrinos {ys) from a GRB provides 
direct evidence for the acceleration of hadrons into the EeV ra nge, and of photo-pion interac- 
tions in the GRB. The GRB standard model (jMeszarosll2006l . and references therein), which 
is based on internal/external shock acceleration, has been used to describe the general fea- 
tures of a GRB and the observed multi-wavelength afterglow. However, the standard model 
cannot well reproduce recent observational results. The early X-ray afterglow s detected by 



Swift exhibited a canonical behavior of steep-flat-steep in their li ght curve (INousek et al. 
20061 ). In 10%-15% of GRBs, precursor activities were observed (IBurlon et al.l 120081 ) . In 



some cases, the precursor preceded the main burst by several hundred seconds with signifi- 
cant energy emission. In the Fermi observa tions of GRB0905 10 and GRB090902B, spectral 
fits revealed a hard power-law component ( lAbdo et al.l 120091 ). The Swift and Fermi obser- 
vations of GRB0905 10 detected gamma rays (7s) in the GeV range up to 200 s after the 
lower energy trigger ( iDe Pasquale et al.ll2010l ). Although many authors have proposed the- 
oretical models to reproduce the complicated time evolution of GRBs an d the high energy 



comp onents in the prompt emission, none of these models are conclusive ( lAckermann et al. 



20101 ). To better understand the ambiguous mechanisms of GRBs, observational probes of 
the optically thick region of the electromagnetic components, as well as hadron acceleration 
processes throughout the precursor, prompt, and afterglow phases are required. Very high 
energy (VHE) vs can be used as direct observational probes, which are effective even in op- 
tically thick regions. A monitor search with sufficient time and spatial resolution and survey 
capability for VHEz/s associated with GRBs is plausible. 

The earth -skimming tau neutrino (z/^) technique, which detects extensive air showers 
( lFargionll2002l ). has the advantage of a large target mass, since it uses air showers produced 
by decay particles of tau leptons (rs) in the atmosphere as the observed signals, rs emerge 
out of the side of the mountain or the ground facing the detector; they are the product 
of interactions between VHE Vr and the earth matter they traverse. Above 1 EeV, air 
f luorescence observati ons based on the earth-skimming z/^ technique have been reported 
( lAbraham et al.ll2008[ ). No air Cherenkov observation was made to date based on the earth- 
skimming Vj. technique with air showers induced by r decays (hereafter referred to as the 
Cherenkov r shower method). But it can achieve sufficient detection sensitivity in the 
PeV-EeV region to be useful in the search for vs originating from hadrons accelerated to 
EeV at astronomical objects. Additional advantages of the Cherenkov r shower method 
are its perfect shielding of cosmic-ray secondary particles, highly precise arrival direction 
determination for primary z/^ and negligible background contamination by atmospheric us 
in the PeV-EeV energy range. 
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2. Ashra Experiment and Observation 



The all-sky survey high-resolution air-shower detector (Ashra) is a complex of unconven- 
tional optical collectors that image VHE air showers in a 42° diameter field o f view (FOV) 



covering 77% of t he entire night sky with a resolution of a few arcminutes (ISasakil 12008 



Aita et al. 


2008a: 


Sasaki et al. 


2008) 



was constructed on Mauna Loa at 3300 m above sea level on Hawaii Island, and includes 



an observatory. Ashra-1 uses electrostatic lenses (lAsaoka fc Sasakill201ll ) in addition to an 
optical system to generate convergent beams, enabling a very low cost and high performance 
image sensor, providing a high resolution over a wide FOV. The electron optics use an image 
pipeline to transport the image from the focal sphere of the refiective mirror optical sys- 
tem. After the light from the image is split, it is transported to both a trigger device and 
high-gain, high-resolution complementary metal-oxide semiconductor image sensor. 

One of the Ashra light collectors built on Mauna Loa has the geometrical advantages of 
not only facing Mauna Kea, allowing it to encompass the large target mass of Mauna Kea in 
the observational FOV, but has also an appropriate distance of ~30 km from Mauna Kea, 
yielding good observational efficiency when imaging air-shower Cherenkov lights which are 
directional with respect to the air-shower axis. Using the advanced features, we performed 
commissioning search for Cherenkov r showers for 197.1 hr between October and December 
of 2008. We served limited 62 channels of photomultiplier tubes (PMTs) as trigger sensors 
prepared for the commissioning runs to cover the view of the surface area of Mauna Kea, 
maximizing the trigger efficiency for Cherenkov r showers from Monte Carlo (MC) study, 
as shown in Fig. [TJ Adjacent-two logic was adopted to trigger the fine imaging, by judging 
discriminated waveform signals from each pixel of the multi-P MT trigger sensor. During 
the search period, ~2 hr before the trigger of GRB081203A (jParsons et al.ll2008l ). GRB 
counterpart (R.A. 15:32:07.58, decl. +63:31:14.9) passed behind Mauna Kea, as viewed 
from the Ashra-1 observatory. Using the same light collector but continuously taking fine 
images of split lights just before the trigger and readout sensors in the image pipehne, we 
set a limit on the light curve as a function of time in the region of 300 s bracketing the 
Swift trigger time for GRB081 203A, including the precursor and prompt afterglow optical 
counterpart ( lAita et al.ll2008b( ). 



3. Analysis 

To investigate the features, selection criteria, detection efficiency, and background rate 
for the observation of Cherenkov r shower images, we generated MC events with pri- 
mary energies of 1 PeV to 100 EeV by 0.5 decade steps, which entered into the rock 
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Fig. 1. — Boundary (large red circle) between the inside (open circle) and outside (hatched 
area) of the FOV of the Ashra-1 light collector, which faces Mauna Kea, and the layout of 
trigger pixel FOVs (blue boxes) for Cherenkov r shower observation. Repositioned array of 
the trigger pixel FOVs (upper blue boxes) to check the detection sensitivity with ordinary 
cosmic-ray air showers at a higher elevation. Firing trigger pixels (thick blue boxes) of an 
observed image of a cosmic-ray air shower readout along the trigger (points). An extended 
portion of the trajectory of GRB081203A counterpart (circular arc), the segment of this 
trajectory used in the v^- search (thick circular arc), the ridge lines of Mauna Kea (red) and 
Mauna Loa (green) mountains, the horizon, and Mauna Kea access road are shown. 



of Mauna Kea uniformly and isotropicall y from a sufficiently large aperture. We used a 
geodetic database around Hawaii island (IMooney et al.lll998l ). and surveyed for ourselves 
the position of the observatory and the terrain of the mountain and its surroundings. To 
study the generation and prop agation of rs in the earth and in the mountain, we used 
PYTHIA (ISjostrand et al.l 120011) to sinaulate the charged current interaction of z/t-s with nu- 
cleons and GEANT4 (lAgostinelli et al.ll2003[ ) for the energy loss due to pair production and 
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bremsstrahlung in r pro pagation. Photonuclear inter action was estimated using the differen- 



1997). TAUOLA 



Heck et al.lll998r ) 



tial crosssection giv en in liver Dutta et al.l (1200 ll ) and lAbramowicz fc Levy! ( 
(jjadach et al.lll993[ ) was used for r decays and CORSIKA (thinning=10~^; 
for air showers induced by r decays. To simulate the Ashra detector, we took into account 
the geometry, hght collection area, mirror reflection, corrector lens transmittance, and the 
quantum efficiency of the photoelectric tubes, so that the fine images corresponding to the 
trigger judgement were simulated in an event-by-event manner. A^^-, the deflection angle 
of r with respect to the primary z/^-, was estimated to be significantly less than 1 arcmin at 
energies above 1 PeV due to the physical processes occurring in the rock of the earth and of 
the mountain (Y. Asaoka & M. Sasaki 2011, to be published). In addition, the reconstructed 
T shower axis can point toward the Ur object within an accuracy of 0?1 if the resolution 
of the image is sufficiently high. Therefore, the Cherenkov r shower induced by PeV-EeV 
Uj. is a fine probe into VHE hadron accelerators once an image is obtained with sufficient 
resolution. Such images can be obtained with the Ashra detector. 

The photometric and trigger sensitivity c alibration of the Ashra light collector was based 
on a very stable YAP(YA103:Ce)-light pulser (IKachanov et al.lll992l ) which was placed at the 
center of the input window of the photoelectric lens image tube mounted on the focal sphere 
of the optical system and illuminated it. Non- uniformity in the detector gain due to the input 
light position was relatively corrected by mounting a spherical plate uniformly covered with 
luminous paint on the input window. To correct for the time variation of the photometric and 
trigger sensitivity because of variations in atmospheric optical thickness, which were mainly 
due to clouds and hazes during the observation period, we performed careful cross-calibration 
to compare the instrumental photoelectric response with the photometry of standard stars 
such as BD-I-75D325 of B-magnitude 9.2, for which the detected images passed through the 
same optical and photoelectric instruments except for the final trigger-controlled readout 
device. We estimated the systematic uncertainty on the basis of our understanding of the 
detector sensitivity to be 30% after applying a combination of the above three complementary 
calibration procedures. 

To validate the detection sensitivity and gain calibration for the Cherenkov r shower, 
we detected and analyzed 140 events of ordinary cosmic-ray air-shower Cherenkov images 
for a total of 44.4 hr in 2008 December using the same instruments used in the Ashra light 
collector, but after rearranging the trigger pixel layout to view the sky field above Mauna 
Kea (Fig. [1]). In the cosmic-ray observation, the trigger pixel layout is centered at zenith 
angle of ~65°. Due to the directionality of the air-shower Cherenkov lights, the photomet- 
ric detection of the lights was strongly dependent upon the impact parameter Rp in 
addition to the primary energy E. From a detailed study of MC simulated proton shower 
events generated with CORSIKA, we obtained a correction function to estimate the ob- 



- 6 - 



served primary energy E as a function of A^^ and the long axis length L calculated using 
Hillas analysis (IHillasI Il985l ). where L was used as an estimator of Rp. As a result, the 
total MC reconstructed energy resolution was estimated to be 62% by evaluating the rms 
of the AE/E = {E — E)/E distribution, for which the error was dominated by the am- 
biguity in _Rp. The same reconstruction procedure was applied both to the observed data 
and the MC data. The observed and MC cosmic-ray flux spectra are shown in Fig. [21 in 
which the MC prediction used the typically observed cosmic-ray flux with E >3 PeV given 
by dN/dE = (3.7 ± 1.1) x 10^.E'~^'°(ni^sr s G eV)-^ fitting to th e cos mic-ray flux observa- 
tions in the knee region (lAntoni et al.l l2005l : lAmenomori et al.l |2008| ) with the power-law 
indices —2.7 {E <3 PeV) and —3.0 {E >3 PeV) assuming continuity at the knee point. We 
quote systematic uncertainties of 30%, mainly due to cosmic-ray flux spectrum observations 
around the knee (20%-28%), and partly due to their parameterization (~14%). Since the 
primary cosmic-ray components are observationally undefined, we present the MC predic- 
tion of cosmic-ray flux spectra, assuming either only protons or irons as the primary cosmic 
rays in Fig. [2l In both cases, the observed data and the MC prediction agreed well on the 
normalization and the shape of the distribution within the expected errors. The estimation 
of the detection sensitivity of the Ashra light collector and the validity of the reconstruction 
procedure were well demonstrated. 

For the Ur search, we used image data acquired using the trigger for 197.1 hr in only 
case of the data st atus defined as good out of the total observation time of 215.8 hr. Using a 
Source Extractor (IBertin fc Arnoutslll996l ). we extracted "clusters" by requiring that of fired 



pixels (Npx) assigned to the cluster satisfied Np^ >4, where accumulated charge (Qpx) in the 
fired pixel (px) was required to be Qpx >6 ADC units. The primary cluster C, defined as the 
cluster with the largest accumulated charge of fired pixels in the cluster (Sc), was chosen as a 
pre-selected shower candidate for the event. To ensure that C was within the effective fiducial 
FOV, the coordinates of the center of gravity (the center of the FOV was (0,0)) {Xc,Yc) 
of C were required to satisfy |Xc>| < 13?8, and {Xc,Yc) could not overlap any obstacles 
such as rocks or the Mauna Loa surface within the FOV. Sometimes, "car light events" were 
triggered when the light beams of a car on the Mauna Kea access road entered the FOV. 
We determined the coordinates of the original position (X^ar, ^car) by averaging {Xc, Yc) for 
events that were apparently identified as car light events in the FOV. The geometric distance 
Dear in the FOV between {Xc, Yc) and (Xcar, ^car) had to satisfy Dear > 0?2. Remaining 



night sky backgrounds (NSBs) were rejected by requiring that Sc > 10 Sc-/, where C 
was defined as the second largest cluster. The NSB rejection was also fairly effective at 
removing contamination from cosmic-ray /i events, in which Cherenkov lights from cosmic- 
ray /is passing through the input window glass (8 mm thick) of the photoelectric lens image 
tube tended to separate into two or more clusters in the same image frame. In the MC study. 
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Observed Energy [PeV] 

Fig. 2. — Observed cosmic-ray flux spectrum (filled box) with bars indicating statistical 
and systematic errors and the MC predictions for proton primary (hatched band) and iron 
primary (shaded band) assumptions. The width of the bands shows the evaluated systematic 
error of 30% of the MC prediction (see the text). 

the total selection efficiency for Cherenkov r showers after satisfying all requirements was 
estimated to be 99.9%, and the expected number of events of cosmic-ray /i contamination 
was 5.8x10"^ events. We evaluated a small residual contamination of the final candidate 
samples from secondary particles in large angle cosmic-ray air showers of 1.3x10"^ events 
using CORSIKA with the curved earth option. From the 197.1 hr of observation, five clusters 
were pre-selected as shower candidates after the fiducial FOV cut. Out of these five clusters, 
three were removed by the car light cut and two were removed by the NSB cut, yielding a 
null result. 



4. Results 

Fig. |3] shows the effective areas for Cherenkov r showers induced by from the 
GRB081203A counterpart with changing energy (E^^) as obtained from an MC study 
that assumed seven original positions on the GRB counterpart trajectory. On the basis of 
the above null result and the estimated effective areas, we placed 90% confidence level (CL) 
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Fig. 3. — Effective area for Cherenkov r showers induced by z/^-s from tlie GRB081203A coun- 
terpart as a function of Ur energy (Eu^) assuming seven original positions on the counterpart 
trajectory (at elevation angles of — 3?62 ~ 0?51) with the Ashra light collector. 

upper limits on the z/^ fluence of precursor and afterglow emissions in the PeV-EeV region, 
for two 3780 s periods, the first between 2.83 and 1.78 hr before GRB081203A and the second 
between 21.2 and 22.2 hr after it, as shown in Fig. HI The same fluence limits are set based 
on observations on the dates of 081122 (2430 s), 081124, 081125 (2000 s), 081126-081205, 
and 081206 (1840 s) for 3780 s each night unless otherwise noted. As sources of system- 
atic uncertainty in the MC estimates of the effective detection areas, we considered the Uj. 
charged current interaction cross section, the energy loss by r in the earth and the moun- 
tain, the geological model around the observation site, and the gain calibration of the light 
collector, which were evaluated to be 50%, 50%, 10%, and 30%, respectively (Y. Asaoka & 
M. Sasaki 2011, to be published). The conservatively combined systematic sensitivity error 
was obtained from the square sum of the above uncertainties as 77%, which affects the result 
of the 90% CL upper limit shown in Fig. HI where we assumed a typical E^^"^ flux to ensure 
unbiased constraints on observationally undefined physical mechanisms of a GRB. 



There are several models which predict VHEz/ flux d uring our observation periods. A 
precu rsor model based on existing supernova before a GRB (IRazzaque et al. I l2004ilvietri fc Stellal 
19981 ) allows VHEi/ emission several hours before a GRB from a choked jet generated by 
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Fig. 4. — Ashra 90% CL limit (thick continuous hne (red)) on 1/^ fluence of precursor and 
afterglow emissions from the GRB 081203A counterpar t in the PeV-EeV region (see the text). 
For comparison, IceCube limits ( lAbbasi et al.l 120101 ) in the prompt (long d ashed dotted 



(gray )) and precursor (long dashed dotted (blue)) phases and RICE limits ( iBesson et al. 



20071 ) in the prompt (dashed (green)) and afterglow (dashed dotted (magenta)) phases are 
shown. The energy ranges were defined as containing 90%, 90%, and 80% of the expected 
signals from assumed spectrum for Ashra, IceCube, and RICE, respectively. 



mass accretio n. A cannonball (CB) m odel also predicts high-energy neutrino emission from 
a choke d CB (IDar fc De Ruiulall2001). Af t erglow rnodels based on externa l forward shock in 



a GRB (Vietri 



1998a 



b; Vietri et al. 



20031 iLi et al.l 120021 ) and CB model flDar fc De Ruiula 



20081 : iKringd l2010l ) predict long-lived VHEi/ emissions for more than a day after a GRB. 
To evaluate quantitatively the afterglow flux in external forward shock at 21.2-22.2 hr 
after the GRB, the VH Eu flux expected fro m an afterglow model assuming interstellar 
medium in GRB030329 jRazzaque et al.ll2004h was rescaled to the flux from GRB081203A 



by using isotropic gamma-ray luminosity {L]^°), redshift (z), and luminosity distance {cIl). 
L]^° = 1.3 X 10^^ erg s'^ was calculated as an average luminosity during time p eriod con- 
taining 50% of the gamma-ray emiss ion (IKann et al.ll2010l : lUkwatta et al.l 120081 ). z = 2.05 
was obtained from iKuin et al.l ( l2009l ). Taking into acc ount the expectation that VHEi/ flux 
is proportional to optical afterglow flux (ILi et al.ll2002l ). the VHEz/ flux was corrected by the 
ratio of magn itudes (0.058) mea sured after 1 . 5 hr from each GRB between GRB081203A 
(jVolkovl |2008|) and GRB030329 (jPrince et al.ll2003l ). The use of apparent magnitude cor- 
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rected the effect of cIl difference. In addition, to account for decay of the optical hght curve, 
the VHEi/ flux was further corrected by the ratio of rn agnitudes ( 4.8 x 10~^) betwe en 21.7 hr 
jAndreev et ahlboosl : IlIu et al.ll2008l : iMori et al.l[2008l ) and 294 s JWest et al.lboosh after the 
GRB081203A. Finally, the redshift of ne utrino energy was con sidered to obtain the neutrino 
spectrum. Although the original model (IRazzaque et al.ll2004j ) calculated the VHEz/ flux in 
reverse shock, it was applied to forward-shock afterglow by taking into account the optical- 
afterglow correction. As a result, our limit between 21.2 and 22.2 hr after GRB081203A 
rejects a flux of 2 x 10^ times larger than the model assuming neutrino oscillation. If wind 
envi ronmerit is a. ssumed and the reverse shock model is applied with the co rrection due t o 
dL (IWrightl I2OO6I ) and dependence of the expected neutrino luminosity (IVietril Il998bl ). 
the same limit rejects a flux of 3 x 10^ times larger than the model. Although our result is 
in agreement with the model, it might constrain some extreme cases such as very effective 
proton acceleration and/or pion production in forward shock. Due to the lack of detailed 
knowledge of the evolution of the accelerated proton spectrum in t he forward shocks, there 
is huge ambiguity in calculation of neutrino fluxes (ILi et al.l |2002| ). It is important to set 
observational limit especially on each burst which would have a much different environment 
burst by burst. 

For comparison, Fig. H ] shows other obser vational limits on the u fiuence from point 
sources (lAbbasi et al.ll2010l : iBesson et al.ll2007| ). Our results are the most stringent in the 
PeV-EeV region and complementary to the IceCube results for the sub-PeV energy region, 
and indicate the advanced instantaneous sensitivity of the system even during this commis- 
sioning phase. Full-scale Ashra observations are expected to have 100 times better sensitivity, 
since the trigger pixel size is halved (1/4 the pixel area) and the trigger sensor will cover the 
entire FOV of the light collector. With this higher sensitivity, contamination of 0.55 clusters 
of air-shower secondary particles is expected in one year of observation data. The advanced 
angular determination of Ashra for Cherenkov r showers to within 0?1 will allow perfect re- 
jection against contamination from air-shower secondary particles, and will provide a viable 
search method for earth-skimming z/^ events, fully utilizing the zero background conditions. 
Our first search for PeV-EeV z/^ reported in this Letter complements other experiments 
i n energy ran ge and methodology, and suggests the prologue of "multi-particle astronomy" 
( lSasakill2000[ ) with a precise determination of time and location. 
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the Coordination Fund for Promoting Science and Technology and by a Grant-in-Aid for 
Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology 
in Japan. 
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